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We show that a finite magnetic flux threading a multichannel semiconductor ring induces spin accumulation 
at the borders of the sample when the Rashba spin-orbit interaction is taken into account even in the absence of 
external electric fields. 
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In recent years semiconductor devices in which the spin or- 
bit (SO) interaction plays a significant role, have been pro- 
posed for applications in spin controlled transport 1 . Among 
them, multiple connected mesoscopic geometries are natu- 
ral candidates to explore how spin dependent effects manifest 
in the quantum interference patterns that appear in transport 
measurements at low temperatures 2 ' 3 . In the case of quasi- ID 
rings, the effect of the SO interaction has been addressed theo- 
retically in a series of paper aVi 6 . As a consequence of the SO, 
the wave function acquires a non trivial spin dependent topo- 
logical phase 7 that manifests in several remarkable quantum 
phenomena. 

We will be considering SO coupling of the Rashba type 
which arises on a 2DEG of a semiconductor heterostructure 
due to the inversion asymmetry of the confining potential 8 . 
When a ring is pierced by a magnetic flux, the SO modifies 
the magnetic flux dependence of the spectrum and therefore 
the conductance 9 and the persistent current (PC), in the case 
of an isolated ring, change as compared to the case without 
SOiS. So far the theoretical analysis has been restricted to ID 
geometries and quasi ID in the two band approximation 11 , 
where the results are qualitative the same as in the ID sys- 
tems. The multichannel nature of realistic rings employed in 
the experiments 1213 and recent controversies around the pos- 
sibility of experimentally observe a spin dependent phase in 
transport experiments with many propagating mode a 14 i 15 chal- 
lenges us to address exactly the 2D geometry. 

In this work we show that when a multichannel ring with 
Rashba SO coupling is pierced by a magnetic flux a spin ac- 
cumulation effect is developed on the boundaries of the sam- 
ple. In addition, even for an even number of electrons, a finite 
spin polarization in the direction perpendicular to the plane of 
the ring is generated whose intensity can be controlled with 
the magnetic flux. These phenomena share some analogies 
with the intrinsic Spin Hall Effect (SHE) studied in bar or T- 
like geometries 161718 , but in the present case the system is in 
equilibrium without external electric fields or voltage drops 
applied^. 

We start by considering a 2D electron gas in the xy plane 
confined to a mesoscopic annular region (multichannel ring) 
threaded by a magnetic flux <£>. The single particle Hamilto- 
nian describing an electron of effective mass m subject to the 
Rashba SO coupling reads 
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where a is the strength of the Rashba spin orbit (RSO) cou- 
pling and the Pauli matrices a are defined as standard. Em- 
ploying polar coordinates p and ip the hard wall confining po- 
tential defining the ring is 



V(p) = 



for r < p < R 
oo otherwise . 
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where r and R are the internal and external radii of the ring. 
The vector potential which is introduced in the Hamiltonian 
via the substitution, p = fik — —ihV — -A, is written in the 
axial gauge as A = (Q/2np)(p . Using , a p — cosipcr x + 
sin ipa y and <r v = — sin ipa x + cos (pa y we can rewrite the 
Hamiltonian as 
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where v = ^- is the magnetic flux in units of the flux quantum 

$o = hc/e. As J z = l z + s z = —ihdcf, + \ha z , commutes 
with H, the eigenfunctions can be chosen as 
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where J z ipj — hjipj and j = I + |. In what follows it will be 
useful to work with dimensionless variables. With that pur- 
pose we define the dimensionless coordinate £ = p/R, the 
aspect ratio A = R/r and 
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with the boundary conditions 

^(A- 1 ) - ffi+xCA" 1 ) = = = m+i{l) ■ (6) 

We can look solutions of the form /;(£) ~ Yi- V (k£) and 
~ Yi + i_ v (k£) where Y;(£) are Bessel functions of 
the type </;(£) or Ni(£). The Rashba term simply acts as ris- 
ing or lowering operator on the Bessel function basis since the 
following standard recurrence relations hold 20 : 
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This is indeed the property which allows to obtain, as in the 
case of a disk geometry 21 , an exact analytical solution. Due to 
the RSO, the bulk spectrum has two branches 



e = k 2 ± (3k 
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Therefore for a given value of e there are two non-trivial solu- 
tions for the momentum k that we denote k + and k~ respec- 
tively. It is then possible to obtain a solution as 
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with 
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and with Y3 and Y4 obtained from Y\ and Y2 by exchanging 
Bessel functions of type J by Bessel functions of type N. 
Defining Y as, 
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the boundary conditions lead to the equation det(Y) = 0. 
Given (3, A and v we solve this equation to obtain the (dimen- 
sionless) energies ej^u) where j is the total angular momen- 
tum and i labels the different eigenstates for a fixed j, in such 
a way that for v = and (3 = Owe have e 



< e 



In order to fix numerical estimates for the parameters we 
consider characteristic values extracted from recent experi- 
ments performed on semiconductor heterostructures defined 
on a 2DEG. Rings with external radius R ~ 400 — SOOnm 
and an aspect ratio A ~ 2 have been recently employed as 
devices 13 . Typical values for the Fermi wavelength are Xp ~ 
40 — 50nm that give kp ~ O.lnm -1 . For R = 400nm one 
gets a maximum value of the (dimensionless) Fermi energy 
€p = (Rkp) 2 ~ 1600. For an effective mass m ~ 0.042m e , 
a Rashba coupling constant a = 8meVnm and R ~ 400 we 
obtain /3 = 2000 nm/R(nm) ~ 4. These parameters de- 
fine the sample S studied in the present work. As the relevant 
situation for an isolated system, we work in the canonical en- 
semble keeping fixed the total number of electrons N as the 
magnetic flux is varied. One can estimate N at zero flux 22 , 
that in this case gives N = [3 ep/8] — 600 (the symbol [..] 
denotes integer part). The maximum number of transverse 
channels M can be then calculated as~ 
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Thus for N ~ 200 — 600 one gets M ~ 4 — 8, in agreement 
with the reported experimental values 12 . 

For v = and finite (3, the SO breaks the degeneracy be- 
tween states differing in one unit of j. The degeneracy be- 
tween states with opposite values of j is broken by the pres- 
ence of a finite magnetic flux v, being the charge PC the sig- 
nature of this broken symmetry. 



As discussed previously in the literature, in ID rings the 
RSO induces a topological phase 9 , A# = ?L2± x <r) • til, 
that once added to the Bohm-Aharonov one, 2ttv, leads to an 
"effective flux" 



'eff 



l(l Ty /(2a ma/h 2 ) + 1) (the ± 
sign depends on the sign of the z spin projection in the local 
spin frame and a is the radius of the ring). It is then via this 
effective flux that the SO interaction affects the behavior of 
the PCi£. 

The situation in the 2-D ring is considerably more involved. 
We display in Fig[0 some regions of the spectrum for both 
(3 = and (3 = 4 cases. The upper panel of Fig^shows the 
lowest eigenvalues of the multichannel ring S as a function of 
v. As the first transverse channel is active in that region, the 
spectrum is similar to that of the ID ring (due to the symmetry 
respect to v = 0.5 the spectrum is shown for < v < 0.5). 
We can observe the evolution with v of the single particle en- 
ergy states labeled by the quantum numbers For (3 — 
(bold dotted lines) we have the double degenerate fundamen- 
tal states, (±1/2,0), and the doublets (1/2, 1), (3/2,0) and 
(-1/2, 1), (-3/2, 0). Notice that at v = 0, except the funda- 
mental state, the others are four-fold degenerate and, as in the 
ID regime, crossings occur only at v — 0, 0.5. The effect of a 
finite RSO coupling is clearly visible in the figure. The RSO 
interaction lowers the energy of each state and it can even 
change the order in which they appear. For the case shown 
in FigQ the lowest energy states are (3/2,0) and (—3/2,0) 
that for finite flux remain almost degenerate and appear in the 
figure as a single line. We then have four states which for very 
small v (i.e, before any v induced level crossing) are ordered 
as (-1/2,0), (1/2,0), (5/2,0), (-5/2,0). Higher in energy, 
we display the states (1/2, 1), (-1/2, 1), (-7/2, 0), (7/2, 0). 
Notice that as a result of the RSO new crossings appear. As 
an example, we draw arrows in the panel as guides for the lo- 
cation of the new crossings between (1/2, 0) and (5/2, 0) and 
(—1/2, 1), (—7/2,0). These crossings are indeed the finger- 
prints of the effective flux v e f f mentioned above^i. 

When many transverse channels are activated the spectrum 
displays additional crossings between levels belonging to dif- 
ferent channels, even in the absence of RSO coupling (see 
lower panel of Fig[0. The crossings that arise due to the 
SO interaction are mixed with the crossings between levels 
with different transverse channel number and it is not straight- 
forward to identify the signature of an effective flux like in 
the ID case, when only a single transverse channel is ac- 
tive. This can be understood looking at the functional form 
of the 2D Hamiltonian Eq.(|3j, whose last term contains the 
ratio a/p between the RSO constant and the radial coordi- 
nate. Therefore, loosely speaking, on average each transverse 
channel feels a topological phase that depends on the value of 
the transverse quantum number. This argument could be ex- 
tended to explain the difference in patterns of conductance os- 
cillations of single-channel and multichannel open rings with 
RSO interaction 15 . In terms of the dimensionless variables, 
the only non vanishing component of the charge current den- 
sity for eigenstates as given in Eq.@ reads, 
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FIG. 1 : Dimensionless energies e as a function of the magnetic flux 
v for the annular cavity S defined in the text. Upper panel: lowest 
eigenvalues in the absence of RSO (bold dotted lines) and for /3 = 4 
(small dotted lines). The arrows are guides for the eyes to locate the 
position of the crossings due to the RSO coupling. See the text for 
details. Lower panel: high energy region showing many crossings 
between different transverse channels, in the absence of RSO (bold 
dotted lines) and for f3 = 4 (small dotted lines). 




FIG. 2: Persistent current 7 as a function of v in the absence of SO 
interaction (dashed lines) and for /3 = 4 (solid line). Upper panel 
corresponds to N = 6 and lower panel to N = 182. The finite slope 
of the jumps is due to size of the flux step (5.10~ 3 ) employed in the 
numerical calculations 



Employing the probability and spin densities, 6j (£) = ff + 
^a z ^ and (<r p )j = ^& p ^, the current can 



be written as 

eh 

v = ^m, 

(13) 

Therefore the effect of the SO interaction in the charge current 
density is unveiled in the last term of Eq. fot . To calculate the 
total charge PC we have to sum the contributions of all states 
up to the Fermi energy, 



(14) 



ometrical factor that takes into account the area of the outer 
circle of the sample I is the magnetic moment. In terms of the 
dimensionless variables, 



}= cEo I= _d± 
4> dv 



(15) 



For the parameters quoted before for sample S, results I ~ 
0.2 I(nA). In Fig|2]we plot I as a function of v for N = 
6(M = 1) and N = 182 (M = 4) in order to show how 
active open channels modify the behavior of the PC and the 
determination of an effective flux as in the ID regime 25 . 

As a result of the RSO interaction, spin projection is not a 
good quantum number. For N occupied states, the mean value 
of the z-projection of the spin is proportional to 
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where as before, n labels the occupied states. Besides a ge- 



where n — 1, N(cf) labels the occupied states and er z (£) cor- 
responds to the total spin density for N particles. Even in the 
case of N even, a non zero spin density is obtained when SO 
is present and v ^ 0. In addition a spin accumulation ef- 
fect is developed and it is manifested in the tendency of a z to 
be positive on one border of the sample and negative on the 
other one. Although the accumulation becomes stronger as 
N and the number of open channels increase, it is also present 
when only one transverse channel is active. In order to explain 
the origin of the effect we first concentrate in a pair of eigen- 
states that, being degenerate at v = with opposite value 
of j and with opposite out of plane spin projection, become 
non-degenerated for finite v. The expectation value of a z in a 
given is (a z )j = \fi-v\ 2 - \gi- v +\ ?■ It is straightforward 
to verify that at v = states with opposite value of j have op- 
posite value of (er 2 ) and therefore for even number of particles 
TV = 2p , is <7 Z (£) = 0. As v > 0, two facts induce a spin ac- 
cumulation effect. On one hand, the magnetic field breaks the 
symmetry between single particle states with opposite value 
of j and on the other hand, due to the presence of additional 
level crossings the total J z of the ground state can change and 
be different from zero for finite flux. In Fig. [5] we show a z 
for p = 4, v = 0.2 and N = 2,4, 6. For N = 2, the single 
particle states have j = ±3/2. As v is turned on, the effective 
orbital index l e ff = I — v becomes 1 — v and —2 — v for 
j = 3/2 and j = —3/2 respectively. Thus as the modulus of 
the effective orbital index decreases (increases) for j = 3/2 
{j = —3/2), the probability density and < a z >j<-j) are 
pushed toward the internal (external) boundary of the sample. 
This symmetry breaking explains the observed accumulation 
effect in this case. For iV = 4 the situation is different. Due to 
the level crossing mentioned before, the single particle states 
in the ground state have j — 3/2, —3/2, —1/2, 5/2 and the 
accumulation effect is mainly due to the unbalance between 
the last two states. We show the case N = 6 which again is 
similar to the case N = 2. 

As the particle number is increased to the relevant exper- 
imental values, transverse channels are activated and the de- 
scription of the effect becomes more complicated. Neverthe- 
less our simulations suggest that the accumulations effect is 
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FIG. 3: Dimensionless total spin density a z as a function of £ for 
v = 0.2 for the sample S. N = 2 (solid line), N = 4 (dashed line) 
and iV = 6 (dotted-dashed line). 
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FIG. 4: Dimensionless total spin density a z as a function of £ for 
v = (solid line) , v = 0.1 (dotted line) and v — 0.45 (dashed line) 
for the sample S. a) TV = 181, b) N = 182. Note that in this case is 
a z — for ^ = as expected for an even number of particles. 



a generic feature of the system in the presence of a magnetic 
flux. As an illustration, in Fig. 0]we show the spin density a z 
as a function of the dimensionless coordinate £ for AT = 181 
and N = 182 and different values of the magnetic flux in the 
range < v < 0.5. For N = 181 and v = (see FigH 
a)) a finite a z (£) is obtained which value corresponds to the 
last occupied level, as expected for odd particle number. On 
the other hand, for N = 182 is cr z (£) = at v = 0, as it is 
shown in Fig b). For v > the spin density profile has a 
more complicated structure than in the quasi ID regime due 
to the behavior of the radial components of the spinors as the 
number of transverse channels increases and the flux is varied. 

In summary we have shown that a finite magnetic flux 
threading a multichannel semiconductor ring with SO interac- 
tion induces spin polarization with opposite sign for the two 
borders of the ring. We believe this system constitutes a new 
proposal to detect accumulation effects induce by SO interac- 
tion in constrain geometries in equilibrium, that is without ap- 
plied electric fields or currents. The characteristic wavelength 
of the accumulation effect is of the order of .1/im, which is 
not far from the sensitivity of methods employed recently to 
probe spin polarization in semiconductor channels 16 . Besides 
the accumulation effect, the integrated spin density S 2 is dif- 
ferent from zero and is sensitive to the value of the magnetic 
flux, as can be inferred from Fig|4] In the presence of an exter- 
nal magnetic field perturbation, the spin magnetization should 
be proportional to £„. With the help of new experimental 
techniques based on resonant methods it should be possible to 
sense changes in the total magnetization of isolated rings due 
to SO interaction^ 
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